Unified treatment of hyperfine splittings in bb, cc and ss. 
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We present a relativistic potential model for bb, cc and ss using the Cornell potential to extend a 
two fermion wave equation determined for the study of the hyperfine structure of hydrogen atom and 
positronium. The hyperfine splitting is obtained by the first perturbation order of the appropriate 
Breit term. The excellent agreement of the numerical with the experimental results for such a 
simple potential shows the fundamental importance of relativity in potential quark models. Only 
three parameters are used for each system. Moreover the same value for the string tension parameter 
has been used for fitting both the bb and cc spectra. The numerical results for the ss system show 
the large infiuence of the hyperfine splitting on the spectrum and are in agreement also with some 
of the not assigned states. 
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INTRODUCTION 



Potential models of interacting qq systems have a long 
history and are still a very lively subject of investigation: 
this is witnessed by the large number of research papers 
and reviews that keep being published which we re- 
fer to for bibliography and exhaustive details on the sub- 
ject. For systems composed by heavy quarks the starting 
point is very often the Schrodinger equation with a po- 
tential having a Coulomb behavior at the origin and con- 
fining at infinity, as the well known Cornell potential; the 
relativistic corrections together with the spin-orbit and 
the spin-spin contributions, thought to be determined by 
one gluon exchange, are taken into account by modify- 
ing the potential with the appropriate correcting terms. 
These are usually treated using a perturbation expan- 
sion in the inverse of the heavy quark mass and give the 
splitting for the states of angular momentum multiplets. 
Attempts have also been made to overcome the limita- 
tions of a potential model due to asymptotic freedom at 
short distances and to light quark creation: a descrip- 
tion of these effects has been tried by means of screened 
potential softening the Coulomb interaction at the ori- 
gin and by letting the confining term saturate at infinity. 
The spin dependent interactions are then modeled by the 
Breit-Fermi potential with a 5-function centered at the 
origin, which may yield some difficulties in explaining 
the hyperfine splittings in the spectra. Although this 
approximation may be good for heavy mesons, a smear- 
ing of the (5-function has been proposed in order to get 
a better description of the small distance behavior. Re- 
cent results 0, however, show that this point is not yet 
firmly settled. We should also mention that the need of 
a more complete and consistent relativistic approach to 
the quarkonium spectrum has sometimes been declared 
and relativistic models have been actually worked out 



[3|. The proposed treatments are often connected with 
field theory along the lines of the Bethe-Salpeter equa- 
tion and the spectra of the resulting equations are not of 
straightforward computation. 

In this paper we present an approach to potential mod- 
els complementary to the lines so far sketched. We use 
the simplest Cornell potential, focusing instead on a com- 
pletely covariant formulation and on the fermionic nature 
of the elementary constituents. Due to the similarity be- 
tween one-gluon and the one-photon exchanges, we also 
add a Breit term: it arises from the expansion in the in- 
teraction delay of the spatial part of the vector potential 
and it describes the so called chromomagnetic interac- 
tion. Due to its origin, this term must be treated at a 
first perturbation order. Further improvements of the po- 
tential are an important issue which should be developed 
at a more phenomenological level of the investigation. 

Our formulation originates from a wave equation for 
two relativistic fermions with arbitrary masses This 
was obtained in a simple quantum mechanical scheme 
by coupling two Dirac equations, each one describing a 
charged fermion in the electromagnetic field generated 
by the other. In the first of 4] we studied in detail the 
pure Coulomb problem proving its covariance, the cor- 
rect Schrodinger limit as well as of the one-particle Dirac 
limit when one of the masses tends to infinity; we also 
proved the cyclicity of the relative time that avoids the 
difficulties of the relative energy excitations. The ra- 
dial equations were obtained and the two-body Coulomb 
spectrum was found. In the second paper [J] the hy- 
perfine structure was investigated by adding the vector 
potential contribution. The spectrum of Positronium in- 
cluding the hyperfine splitting was then calculated. We 
also got numerical evidence confirming that the spin-spin 
Breit interaction should be treated at the first perturba- 
tion order: including higher orders or even the complete 
non-perturbative contribution gives back wrong answers, 



as expected. Eventually we compared our results with 
those obtained by QED semi-classical expansions finding 
an agreement better than up to the fourth power of the 
fine structure constant. 



THE TWO-FERMION WAVE EQUATION WITH 
CORNELL POTENTIAL AND BREIT TERM 



Although we cannot reproduce here the derivation of 
the two-fermion relativistic wave equation given in [3] 
we want, however, to stress some of its properties. It 
is conceptually very simple; the calculations to find the 
radial system, although rather long, are straightforward 
and modeled on the scheme of the spherical spinors. The 
Dirac equations entering the wave equation prescribe, 
then, the correct form for the interactions according to 
their tensorial nature: the Coulomb-like term of the Cor- 
nell potential is vectorial and thus minimally coupled to 
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TABLE I. The bb levels in MeV. In the first column the 
term symbol, the I'^{JPC) numbers and the particle name. 
cr=l.lllGeV"^, q:=0.3272, m;, =4725 . 5 MeV. Experimen- 
tal data from (6(1 . 



the energy; the linear term is scalar and therefore coupled 
to the mass. We recall that only a scalar growing po- 
tential is actually confining, while an unbounded vector 
interaction is not [Sj . Hence in our case case it is inappro- 
priate to write the Cornell potential just as the sum of the 
two terms as in the Schrodinger limit where the distinc- 
tion becomes meaningless. The equation for the qq model 
is obtained by using [4] and coupling the scalar interac- 
tion to the fermion masses. We call r^, qa the Wigner 
vectors of spin one given by the spatial parts of relative 
coordinates and momenta boosted in the frame with van- 
ishing total spatial momentum and we put r — {raTaY^"^ 
(sum over repeated indices). The vector and scalar cou- 
plings produce the terms E + a/r and m + {l/2)ar. The 
final form of the two fermions equation contains thus only 
the three parameters a, a, m and reads 

(7(1)7(1), - 7("2)7(2) J^a + (7(°i)+7['2)) ("^ + \^ 

- (s+ ^) + V^s(?-)] = 0. 

In ([T} 7(j) are the gamma matrices acting in the spinor 
space of the i-th fermion and 



(1) 



^B{r) = ^7('l)7(l)„7('2)7(2)b(^a6- 



(2) 



denotes the Breit term generating the hyperfine split- 
ting. We found in jj] a very convenient way to evaluate 
the first order contribution of this term without resort- 
ing to eigenf unctions: substitute in ^ VB(r) by eVsir) 
and calculate the first derivative of the eigenvalues with 
respect to e in e = from the numerical solutions of 
the differential equations. We also proved that the use of 
eigenfunctions would, obviously, produce the same result. 

The radial system is obtained by diagonalizing angular 
momentum and parity. As in [J] it is formed by four 
algebraic plus four first order differential equations for 
each parity. Using the algebraic relations and defining 
the dimensionless variables , w , s by 



E ^m{2 + nw), 



(3) 



the radial system for ([T]), replacing Vsir) by eVB(r), is 
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Here Ao = yle|e=o ^0 = Se|e=o and u' = du/ds. Letting 
J^=j{j + 1), the coefficients for the even parity are: 
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TABLE II. The cc levels in MeV. cr=l . Ill GeV" 
mc = 1 394. 5 Me V. Experimental data from [^. 



■ , a=0.435, 



with/i(s) = {2 + VLw)/\/n + a/s, k{s) = {2 + VLs) / {2^/^^). 
The coefRcients for the odd parity systeni are: 
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Before presenting the resuhs, a word is in order about 
the numerical method we have used. The solution of the 
singular boundary value problem has been obtained by a 
double shooting method, the spectral condition being the 
vanishing of the 4x4 determinant of the matching condi- 
tions in a crossing point Pade techniques have been 
used in order to improve the accuracy of the approximate 
solutions in the neighborhood of the two singular points, 
zero and infinity. The integration precision has always 
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TABLE III. The Breit correction As in MeV for some levels 
of bb, cc, ss. 



been kept very high and tested against the stability of 
the spectral values. 



NUMERICAL RESULTS FOR bb AND cc 

Let us present the numerical results for the levels of 
bb and cc. In order to have a test as good as possible of 
the relevance of the relativistic dynamics in quarkonium 
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TABLE IV. The ss levels in MeV. cr=1.34 GeV" 
ms=134 . 27 MeV. Experimental data from 



a=0.6075. 



3 



models, we have aimed at choosing the least number of 
parameters for fitting the experimental data. Flavor in- 
dependence could be expected for heavy quarks. In fact, 
doing separate fits for bb and cc we find that the string 
tensions turn out to get the same value within the com- 
putation precision: we can thus directly assume a single 
parameter a for both models so that the total number 
of parameters is only five. The parameters a, constant 
in each system, are independent on the energy scale and 
they are determined by separate fits. However the ratio 
a{bb)/a{cc) = 0.752, numerically found, is very close to 
the ratio of the two as{r]b{lS))/as{r]cilS)) ^ 0.746 for 
average values A5 = 0.30, 0.33 GeV for n/ =4,3 ^. 

The spectra show common features, generally shared 
by all potential models: the states group into doublets of 
s states and quadruplets of p, d, ... states. It clearly ap- 
pears that the results are in very good agreement with ex- 
perimental data below the threshold of B and D mesons 
[6| for bb and cc respectively. Above the threshold the 
calculated energy of the levels is larger than the experi- 
mental one and a softening of the potential growth could 
make a sensible difference in the agreement with the data 
for the higher levels. The regularity of the level pattern 
is however maintained. For the resonance X(3782), hav- 
ing the two possible assignments J^*^ = 1++ and 2 ^ 
0, the model could indicate a Xci classification. Noth- 
ing can be suggested for X(3915) and X(3940), having 
no accepted quantum numbers. 

The situation is simpler in Table I, where there are 
no unclassified physical states. We point out the good 
estimate of the recently discovered Xbi^P) resonance 
0, staying just below the B production threshold. In- 
stead, as expected, the calculated values for T(4^si) and 
T(5^si) exceed the experimental values. 

We present in Table III the calculated absolute values 
of the Breit corrections for some states of bb, cc and ss: 
they are responsible for the hyperfine splitting and have 
to be subtracted from the initial pure Cornell potential 
levels in order to obtain the values given in Tables I, II 
and IV. As expected, the corrections decrease for increas- 
ing values of J and become more and more important for 
decreasing masses of the component quarks. 

A LOOK AT ss AND SOME CONCLUSIONS 

The good agreement obtained for bb and cc has urged 
us to test our model also on the more elusive ss system, 
for which there are few accepted experimental states and 
many puzzles that cannot be explained by a simple po- 
tential model. We believe, however, that these numerical 
results may be of interest, since the much lighter mass 
of the s quark highly enhances the relativistic character 
of the ss composite system and the fundamental role of 
the Breit corrections, giving rise to huge hyperfine split- 
tings. Due to these reasons the string tension a has not 



been given the same value of the previous systems but 
has been considered a fitting parameter, finding a value 
larger than in bb. We report our results in Table IV, 
where we have also included the unassigned fi(1420), 
X(1750), (/)3(1850) and 0(2170) . Although we cannot 
have a complete phenomenological confidence in the nu- 
merical results, still a fair number of experimental data 
can be accommodated with a pretty good accuracy. For 
instance the model could suggest a (l^di) assignment for 
X(1750). 

To conclude we can summarize that the two-fermion 
relativistic equation describes in a unified way three dif- 
ferent types of qq systems with the same treatment of 
the hyperfine splitting as the one successfully applied to 
positronium. With a very limited number of parameters 
it provides satisfactory fits for the qq levels below the 
thresholds of lighter meson production and a concrete nu- 
merical method to evaluate the effect of the Breit term. 
It appears thus to be an effective tool in the analysis of 
two fermion bound states in a given range of energy with 
a prescribed interaction law. 
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